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Spectropolarimetry of the Type lib Supernova 2001ig 
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ABSTRACT 

We present spectropolarimetric observations of the Type lib SN 2001ig in 
NGC 7424; conducted with the ESQ VLT FORSl on 2001 Dec 16, 2002 Jan 3 
and 2002 Aug 16 or 13, 31 and 256 days post-explosion. These observations are 
at three different stages of the SN evolution: (1) The hydrogen-rich photospheric 
phase, (2) the Type II to Type lb transitional phase and (3) the nebular phase. 
At each of these stages, the observations show remarkably different polarization 
properties as a function of wavelength. We show that the degree of interstel- 
lar polarization is 0.17%. The low intrinsic polarization (~ 0.2%) at the first 
epoch is consistent with an almost spherical (< 10% deviation from spherical 
symmetry) hydrogen dominated ejecta. Similar to SN 1987A and to Type IIP 
SNe, a sharp increase in the degree of the polarization (~ 1%) is observed when 
the outer hydrogen layer becomes optically thin by day 31; only at this epoch is 
the polarization well described by a "dominant axis." The polarization angle of 
the data shows a rotation through ~ 40° between the first and second epochs, 
indicating that the asymmetries of the first epoch were not directly coupled with 
those observed at the second epoch. For the most polarized lines, we observe 
wavelength-dependent loop structures in addition to the dominant axis on the 
Q — U plane. We show that the polarization properties of Type lib SNe are 
roughly similar to one another, but with significant differences arising due to line 



^Department of Astronomy and McDonald Observatory, The University of Texas, 1 University Station 
C1402, Austin, Texas 78712-0259, U.S.A.; jrm@astro.as.utexas.edu; wheel@astro.as.utexas.edu 

^ESO - European Organisation for Astronomical Research in the Southern Hemisphere, Karl- 
Schwarzschild-Str. 2, 85748 Garching b. Miinchen, Germany; fpatat@eso.org; dbaade@eso.org 

^Department of Physics, Texas A&M University, College Station, Texas 77843-4242, U.S.A.; 
wang@physics.tamu.edu 

''Department of Physics, Florida State University, Tallahassee, Florida 32306-4350, U.S.A.; 
pah@astro . physics .fsu.edu 



- 2 - 



blending effects especially with the high velocities observed for SN 2001ig. This 
suggests that the geometry of SN 2001ig is related to SN 1993J and that these 
events may have arisen from a similar binary progenitor system. 

Subject headings: supernovae: general — supernovae: individual: (SN 2001ig) 
— techniquesipolarimetric 



Introduction 



Asymmetries are thought to be inherent with core-collapse S upernovae (CCSN e) and 
are inferred thro ugh observations of runaway O stars and pulsars (IDray et al.ll2005l ). spec- 
tral line profiles (IMazzali et al.ll2005l ) and, in the extreme, through the lii ik between CCSNe 
and highly-coUimated jets that give rise to Gamma Ray Bursts (GRBs; IWoosley fc Bloom 
20061 ) . Spectropolarimetry has been used to directly measure t he geometries of b oth CC- 
SNe JWang et al.ll2003bl : Iheonard et aPbooeh and Type la SNe JWang et al.ll2003ah . While 
J2005h provide an example of studying asymmetries in CCSNe by exam- 



Mazzali et al 



ining the shape of line profiles in ordinary spectroscopic observations, spectropolarimet- 
ric observations can provide a straightforward and unambiguous measure of asymmetries 
(IShapiro fc Sutherland! Il982l ). The former technique can only be conducted at late times 
during the nebular phase, when the ejecta is optically thin and the spectrum is dominated 
by broad emission lines of fordidden transitions. At this later epoch the profiles of these 
lines can trace the geometry of the ejecta, but are dominated by the Circumstellar Medium 
(CSM)-ejecta interaction and are no longer necessarily an accurate tracer of the geometry of 
the explosion mechanism. In addition, during the e arly, optically-thick phases the profiles of 
lines are only weakly dependent on the geometry (IHowell et al.ll200ll ). Alternatively spec- 
tropolarimetry can be conducted at any epoch, especially at early times when the geometry 
of the ejecta is dictated by the shape of the progenitor and the explosion mechanism. 
For a spherically-symmetric, electron scattering atmosphere, the polarization vectors of the 
light received, from different portions of the SN, cancel out leading to no net observed po- 
larization. For geometries that depart from the spherical symmetry, there is incomplete 
cancella tion of the polarization vectors which leads to a net observed linear intrinsic polar- 
ization Jshapiro fc Sutherland! Il982l : iMcCaUlllQsi iMcCall et al.lEosi ). Electron scattering 
is wavelength independent and the measured level of continuum polarization is associated 



^Based on observations made with ESO Telescopes at the Paranal Observatory, under programmes 68. D- 
0571 and 69.D-0438. 
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with the asymmetry of the photosphere. In addition, the polarization associated with spec- 
tral lines, through line scattering, yields information concerning the geometry of the line 
forming region, and specifically the distribution of particular elements within the SN ejecta 
dWang et al]l2003bh . 

SNe are principally classified by the absence (Type I) or presence (Type II) of hydrogen 
in their early time spectra. All types of SNe except one (Type la) are associated w i th the 
core-collapse induced explosion at the ends of the lives of massive stars. I Wang et al.l (120011 ) 
showed that, in general. Type II SNe are polarized at much higher levels compared with 
Type la SNe; at ~ 1% and ^ 0.2 — 0.3% respectively. Spectropolarime tric studies of the 
pecuhar Type II SN 1987A Jjeffervlll99ll ) and the Type IIP SN 1999em JWang et al.ll2002h 
showed a significant increase, from 0.1% to 1%, in the level of the polarization as the time 
since explosion increased, suggesting increasin g asymmetries at inc reasing depths into the 
SNe. In the case of the Type IIP SN 2004dj (ILeonard et al.l l2006l ) the polarization prop- 
erties were directly correlated to the evolution of the SN light curve. During the optically 
thick phase, the polarization remained at a nearly constant level; increasing dramatically 
when the SN entered the optically thin stage, when the complicated internal structure of 
the core was revealed. Type Ibc SNe demonstrate significantly higher levels of polarization. 
Spectropolarimetric observati ons of the Typ e Ic SN 2002ap showed intrinsic continuum po- 
larization levels > 1%, which [Leonard et al.l (120021 ) suggest arises from an asymmetry of at 
least 20%. In addition, the differences in polarization angles between features due to Fe and 
other elements dern onstrated the different distributions of these elements within the ejecta. 
Wang et al. J2003bh compared the O I 7774A line with F e II lines to show that these elements 



had very different geometries, and [Leonard et al.l (12002! ) similarly studied the difference be- 
tween Ca lines and those of Fe group elements. 

The case of spectropolarimetric observations of Type lib SNe is particularly interesting, as 
this type of SN is the link relating H-rich and H-deficient CCSNe. At early times Type lib 
SNe show hydrogen in their spectra, but these features weaken over a period of weeks, leadiri g 
to a hydrogen deficient Type lb SN (IFilippenko. Matheson fc Hdll993l : ISwartz et al.lll993l ). 
The spectroscopic evolution of Type lib SNe is due to the presence of a thin veil of hydrogen 
surrounding the core of the progenitor, which is observed at early times but is subsequently 
dispersed. The princ ipal example of spectro pola rimetry of a Type lib SN comes from obser- 
vations of SN 1993J. iTrammell et al.l (119931 ) and lTran et al.l (119971 ) observed polarizations of 
>, 1%, requiri ng major to minor axis ratios of the emitting region to be between 1.54 and 2. 



HoflichI (119951 ) modelled these polarizations as being due to asymm etries of the outer H-rich 



envelope and an off-centre energy source, and iTran et al.l (119971 ) presented the case of an 
asymmetric circumstellar medium. This is particularly important for the understanding of 
how t he core-collapse mechanism can give rise to a range of SN types. IWang fc Wheeler 
( 119961 ) also showed how, in the case of SN 1987A, the technique of spectropolarimetry may 
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be used to probe the nature of the dust in the immediate vicinity of SNe at late times. 
The observations of the Type lib SN 2001ig and the reduction of the data are given in Sec- 
tion [21 The resuhs of these observations are presented in §31 and are discussed in §11 in the 
context of SN 200 lig as a CCSN and in the broader context of previous spectropolarimetric 
observations of other CCSNe, in particular of Type lib SNe. 



2. Observations and Data Reduction 



SN 2001ig was discovered by lEvansI (120011 ) on 2001 Dec 10.43UT, in t he galaxv of 
NGC 7424. SN 2001ig is located at ^2000 = 22i^57'^30!7 and ^2000 = -41°02'26" jRvder et al. 



2OOII ). which corresponds to an offset of 185"E and 108"N from the centre of the host galaxy. 



The p osition of SN 200 lig, relative to its host galaxy, is shown as Fig. [H iRyder et al. 
(120041 ) estimated an explosion date of 2001 Dec 3 (JD452247.5), fror n models of th e radio 
light curve; we adopt this as the reference date for our observations. iPhillips et al.l (1200 ll ) 
provisionally ident ified SN 2001ig as being of Type lib, resembling SN 1987K. Furthermore, 
Clocchiattil (|2002| ) observed the spectroscopic transition from the H-rich phase to the H- 
deficient Type l b phase, confirniing the classification of Type lib. Foreground Galactic 
reddening, fro m Schlegel et al. ( 1998 lR is E{B — V) = 0.01 or Ay = 0.03 assuming a 
Cardelli et al.l (119891 ) G alactic reddening l aw. T he heliocentric recessional velocity of the 
host galaxy is given by iKoribalski erall toO^ as 939km s"^ . HyperLEDA| quotes the 
recessional velocity, corrected for infall towards the Virgo cluster, to be fvir = 758km s~^ ; 
for Hq = 75km s~^ this corresponds to a distance of lOMpc. 

Spectropolarimetric observations of the Type lib SN 200 lig were conducted at three epochs: 
2001 Dec 16.1, 2002 Jan 3.1 and 2002 Aug 16.3. These observations were made with the Very 
Large Telescope (VLT) at the European Southern Observatory (ESO), P aranal, Chile; with 



the F ocal Reducer/Low Resolution Spectrograph 1 (FORSl) instrument (lAppenzeller et al. 
19981 ). A log of these observations is presented as Table [H Throughout this study all 



times are a given as UT. At each epoch FORSl was mounted on Melipal (Unit 3) VLT 
at the Cassegrain Focus. FORSl, in the spectropolarimetry mode, functions as a dual 
beam spectropolarimeter. The FORSl instrument was used with the standard resolution 
collimator, providing a plate scale of 0.2" px~^. The standard "striped" PMOS slit mask 
was used, with slit width 1" and length 22", and for each set of observations the slits were kept 
at PA = 0°. FORSl uses a sequence of a "super achromatic" retarder plate and a Wollaston 



^http: / / ncdwww.ipac.caltcch.edu 
■^http:/ /leda. univ-lyonl.fr/ 
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prism (ESO No. 34) to separate the different polarization components. Observations of SN 
2001ig were conducted with the retarder plate at 4 angles: 0?0, 45?0, 22?5 and 67?5. At each 
retarder plate position the WoUaston prism produced two beams that were dispersed using 
the grism G300V, providing a dispersion of 2.6A px~^ and spectral resolution, measured 
from arc lamp calibration exposures, of 12. 3A. The observations covered the wavelength 
range 4450-8600A, when the GG435 order separation filter was used, and 3600-8600A, when 
no filter was used. FORSl uses a 2048 x 2048 Tektronix CCD detector, and the observations 
were conducted with gain 0.71e~ADU~^ and readout noise of 5.6e~. The data were reduced 
in the standard manner using and a series of our own routines. The data were 

corrected for bias and overscan. A master unpolarized normalised fiat was constructed from 
fiat observations conducted with the retarder plate at each of the four angles, at each epoch. 
The fiat was applied to object frames, and object spectra, for the ordinary and extraordinary 
rays at each of the four retarder plate angles, were optimally extracted. These spectra were 
wavelength calibrated through comparison with observations of HeHgArCd arc lamps. The 
Stokes parameters Q and U a nd total polarizat ion p and polarization angle 9 were calculated 



in the standard manner (see iJehin et al.l 120051 ). with the data re-binned to ISA to improve 



the signal-to-noise (S/N) for the individual Stokes parameters. Here we have used a flux 
binning technique, with bins recalculated prior to the calculation of the Stokes paranaeters . 
This yields identical results to the weighted Stokes parameter binning of IWang et al.l (119971 ). 
The bin size was determined by the spectral resolution of t he observations, ari d was required 
to be much less than the velocity widths of the SN features (IWang et al.ll2003al ). A correction 
was applied for the wavelength dependent chromatic zero angle offset. The total flux spectra 
of SN 2001ig were flux calibrated using observations of flux standard stars with the retarder 
plate set to 0?0. The data were corrected for the heliocentric recessional velocity of the host 
galaxy. 



3. Observational Results 
3.1. Flux Spectra 

Flux spectra of SN 200 lig, at the three epochs, are presented as Fig. [2J Velocities at 
the absorption minima for different lines, in the spectra of 2001 Dec 16 and 2002 Jan 03, are 
given in Table [2l 



''IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation - ,http: / /iraf.noao.edu/) 
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The spectrum of 2001 Dec 16, at 13 days post-explosion, is dominated by the strong Ha P 
Cygni profile, with an absorption minimum at —19 000 km s~^ . The peculiar profile of the 
line, with sharp edges separating the emission and absorption components, and a redward 
shoulder to the emission component suggests blending with He I 6678A. This interpretation 
is further validated by the presence of an absorption "notch" at 6255A, corresponding to 
He I 6678A also blue shifted by -19 000 km s^^ . There may also be a contribution to the 
Ho; profile due to blended Si II 6355A, which is seen in the spectrum of the next epoch's 
observation, but this is not resolved at the first epoch. An absorption at 5603A is identified 
as He I 5876A (-14 000 km s~^ ), however the associated emission feature is not symmetrical 
about the rest wavelength, the peak is blue ward of the rest wavelength. It may be that the 
emission feature is truncated by some other feature that lies redward. An absorption feature 
at the red end of the spectrum is identified as the absorption associated with the Ca II IR 

triplet ( 12 700 km s"^ ) - with the emission feature redward of the wavelength limit of 

the observation. This absorption is, however, blended with a broad telluric feature, as indi- 
cated on Fig. [21 increasing the uncertainty on the measured velocity at the true absorption 
minimum. There is a single absorption at the rest wavelength of Na ID at the recessional 
velocity of the host galaxy. The measured equivalent width of this feature was determined 
to be Wx(NaID) = 0.7 ± O.lA. Using the relations of iTuratto et all fl2003l ). this equivalent 
width corresponds to two possible reddenings of E{B — V) = 0.10 or 0.32, which are due to 
internal reddening in the host galaxy. 

The spectrum of 2002 Jan 3, 31 days post-explosion, shows the SN at a transitional state 
from a Type II to Type lb SN, with the He lines growing in strength compared to the hy- 
drogen lines. The contamination of Ha by lines from other species is more evident. The 
sharp edge between the absorption and emission components of this profile results from the 
superposition of Si II 6355A , Ha and He I 6678A . As in the previous epoch, the He I 
emission component produces a redward shoulder on the Ha emission component, and leads 
to an absorption "notch" at 6415A (or —12 000 km s"^ ) which suppresses the blue side of 
the Ha emission. The absorption minimum blue ward of the Ha peak is characterised by 
two minima at 6192Aand 6265A , which corresponds to the Si II (—8 000 km s~^ ) and Ha 
(—13 600 km s~^ ), respectively. He I 5876Ais observed as a much stronger feature, with ob- 
vious emission and absorption components; the absorption minimum corresponds to —12 000 
km s~^ . The steep transition blueward from the He I emission to the absorption shows a 
sharp notch due to Na ID. O I 7774A is present as a P Cygni profile with an absorption 
minimum at 7500A, corresponding to —10 600km s~^ . 

The final spectrum, of 2002 Aug 16, 256 days post-explosion, shows the complete transforma- 
tion of the SN spectrum to the nebular phase, dominated by strong forbidd en line emission 



featur es, typical of other CCSNe at the nebular phase (e.g. SN 1993J; iMatheson et al. 



2000b). Most prominent at this epoch is the O I feature at AA6300, 6364A , which shows a 
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number of peaks that are indicative of asymmetries in the expanding ejecta. A crude Gaus- 
sian fit to the O I emission feature gives a full width at half maximum of 6 700km s~^ . Also 
present in the spectrum are features due to Mg I] 457lA, [Ca II] 7291,7324A, and emission 
features due to Fe II 5169A, Hel 5876A and O I 7774A. 



3.2. Spectropolarimetric Properties of SN 2001ig 

3.2.1. General Spectropolarimetric Properties 

The spectropolarimetric observations of SN 2001ig, at the three epochs, are presented 
as Figs. [31 m and O The observed polarization properties are a product of both the intrinsic 
polarization of the SN itself and polarization due to the intervening Interstellar Medium 
(ISM). The calculation of the latter, the Interstellar Polarization (ISP), is an important 
consideration and the determination of its quantity is discussed at length in the following 
sections along with analysis of the intrinsic polarization of the SN. Importantly, the change 
in the polarization properties between 2001 Dec 16 and 2002 Jan 3 is directly correlated with 
the evolution in the spectroscopic features in the SN spectrum. The polarization features 
associated with He I show particularly obvious evolution, with increasing polarization as the 
He I flux line strength increases at the later epoch. In addition, we note that the observed 
polarization P.A. in the first epoch varies widely across ~ 45° in the first epoch, whereas in 
the second epoch the P.A. falls within a narrower range of ~ 20°. The polarization features 
in the second epoch are directly correlated with spectral lines, and for instances with P 
Cygni profiles we se e the classic "inverted P Cygni" profiles in the polarization spectra 



(jCropper et al.lll988l ). These profiles are observed for Ha, He I 6678A and O I 7774A. 
The observations of the last epoch (2002 Aug 16 or 256 days post-explosion) show that, with 
the SN having entered the nebular phase, the spectrum is almost completely unpolarized, 
without significant correlation between polarization features and spectral features. At such 
late epochs it is presumed that geometric dilution of the scattering particles in the SN ejecta 
leads to a drop in the scattering optical depth such that, even in the presence of asymmetries, 
the density of scatterers is not sufficiently high to produce the polarization signature. The 
amount of scatter in the polarization properties of neighbouring wavelength bins is directly 
correlated with the S/N and in instances where the S/N is high, as for the A4571 Mg I], 
AA6300, 6334A [O I] lines and AA7291, 7324 [Ca II] lines, the data has small uncertainties 
and have polarizations tending towards the same values (see §3.2.3p . 
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3.2.2. Data on the Q-U plane 



The Q — U plane is a powerful tool for studying the effects of the ISP and simultane- 
ously examining the behaviours of the polarization and the polarizati on angle with wave- 



length. This tool h as been used to great ef fect by a number of studies (ICropper et al.lll988 



Wang et aliboOSal lbll | Wang et al.l (l2003bl ) showed how the distribution of the data on the 
Q-U plane can be parameterised by a "dominant axis," with functional form: 



U = a + /3Q 



(1) 



This axis can be used to study and quantify symmetries in the polari zation data, which 
are directly linked to physical axial symmetries in the SN. In addition, IWang et al.l (120061 ) 
showed that Q and U Stokes parameters across spectral lines have their own "d ominant" axes. 



leading to information about the relative distribution of particular species. ICropper et al. 



(119881 ) showed that, for spectropolarimetry of Ha of SN 1987A, the progression through the 
points on the Q — U plane, in wavelength order, yielded a loop structure, as opposed to 
a simple straight line. This was interpreted as the symmetry ax is, of the aspheric al SN, 



having different orientations as a function of wavelength (see also iKasen et al 



20031 ). The 



distribution of the spectropolarimetry data of SN 2001ig on the Q — U plane are shown, for 
each of the three observational epochs, as Figs. [6], [7] and [HI 

Since all the observations have been conducted with the slits in the same orientation {PA = 
0°), all the polarization data presented here has Stokes parameters in the reference frame of 
the sky. At the first epoch, the concentration of the data on the Q — U plane is obvious, 
although the location of the dominant axis is not immediately clear. The dominant axis, at 
each epoch, was calculated using a xl minimization technique, that excluded the low S/N 
red end of the spectrum corresponding to the Ca II IR triplet absorption. The low intrinsic 
polarization data at the blue end of the spectrum (< 6OOOA) lies on one side of the dominant 
axis, while the data associated with the Ha feature lie on the opposite side of the dominant 
axis (as shown on Fig. [6]). At the second epoch the dominant axis is clearly defined and has 
rotated through ~ 40° compared to the first epoch. The data is more tightly clustered on 
the Q — U plane, which is partially due to an improvement in the S/N at the second epoch 
as the SN had just reached maximum by day 31. The difference in the locus of the data, on 
the Q — U plane, at the two different epochs results in two different values of f3 (even after 
ISP subtraction, as discussed in Section [3.2.31) . The data at day 256 show no discernible 
structure and are dominated by the statistical uncertainty due to the faintness of the SN, 
compared to earlier epochs, and, hence, low S/N. The use of a dominant axis to describe the 
data at the third epoch was, therefore, inappropriate. 

The dominant axis defines a useful "zero point" for rotating the Stokes Parameters to a 
more meaningful reference frame than that defined by using PA = 0° on the sky. The Q and 
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U param eters can be rotate d using the standard two-dimensional rotation transformation, 
given by (IWang et al.ll2003al ) as: 



Pd = {Q- Qisp) cos 26^ + {U- Ujsp) sin 26^ (2) 

Po = -iQ- Qisp) sin 2ed + iU- Uisp) cos 29 d (3) 

where 9d is the polarization angle of the dominant axis, determined from Eqn. [H and the 
Po is the polarization in the orthogonal direction on the Q — U plane but corresponds to a 
physical rotation of 45°. To ensure rotation about an origin of zero intrinsic polarization of 
the SN, the rotated Stokes parameters are calculated using the Stokes parameters corrected 
for the ISP (see §3.2.31) . The rotated Stokes parameters for the data of 2001 Dec 16 and 
2002 Jan 03 are shown as Figs. [9]and[T0l respectively. 



3.2.3. Determination of the Interstellar Polarization 



The correct study of the intrinsic polarization properties of SN 2001ig requires careful 



accounting for and subtraction of polarization due to the intervening ISM. ISerkowski et al. 
( 119751 ) present the relationship between reddening and the ISP, where Pisp{%) < 9E{B — V) 
for Galactic type reddening and dust polarization laws. This relationship can be used to 
place important maximum limits of the total B/V polarization using reddenings determined 
from non-polarimetric techniques. We assume that dust in the host galaxy N GC 7424 follows 



similar red dening and po l arizat ion laws as the Galaxy, although studies by [Leonard et al. 



( I2OO2I ) and lClayton et al.l (120041 ) suggest that this may not be the case. The maximum polar- 
ization due to the Galactic ISM is 0.09%, given E{B — V) = 0.01. In addition, the reddening 
calculated from the equivalent width of the Na ID lines, in §3.H yields two limits on the 
polarization, due to the ISM of the host galaxy, as 0.9% and 2.9%. Under the assumption 
that the Galactic polarization law holds in NGC 7424, these are strict upper limits on the 
ISP (including the relevant uncertainties in the measurements of the equivalent width and 
in the calculation of the reddening). These limits describe a circle on the Q — U plane, and 
comparison of Figs. [H] to [S] shows that the observed data have polarizations significantly 
lower than these limits. 



Trammell et al.l (119931 ) and lTran et al.l (119971 ) assume that the P Cygni emission components 
of lines such as Ha are inherently depolarized, such that any residual polarization across the 
emission features is due solely to the ISP. At the first epoch. Figs. [3] and Ej, the Ha emission 
component has one of the highest polarization levels across the observed wavelengths. The 
Ha profile, as discussed in Section 13. H is a product of both Ha emission and the blueward 
absorption associated with the He I 6678A. In this case, therefore, the polarization observed 
across Ha is obviously not the ideal case of a completely depolarized emission component of 
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one line and, hence, precludes the use of this assumption for determining the ISP. 
As we have observations at more than one epoch, it is possible to calculate the ISP as a 
non-varying quantity, with a single polarization angle, by comparing the behaviour of SN 
2001ig on the Q — U plane at different epochs. There could also be a local,variable, polarized 
source within the seeing disk of these observations; there is no known technique to correct 
for such phenomena. As noted in §3.2. H the data show a large rotation and an increase 
in polarization between the first and second epochs, corresponding to the change from an 
H-rich spectrum to a transitional He-rich spectrum. The ISP is, by assumption, constant 
between the two ep ochs, requiring our choic e of ISP to be appl icable to both cases. Fur- 
thermore, following iHowell et al.l (I2OOII ) and IWang et al.l (l2003a] ). we expect the ISP to be 
located at one end of the dominant axis of the locus of the data on the Q — U plane and in 
the proximity of the data at bluer wavelengths, due to high degree of blending and, hence, 
depolarization of Fe II lines at bluer wavelengths. For the data of the first epoch (Fig. [H]) we 
note that the Q and U of the lowest total polarization features limit the ISP to Qisp < 0.2, 
Uisp > —0.2. For the second epoch, shown as Fig. [71 the dominant axis is easily identified, 
and the ISP must be located at the low polarization end of the distribution of data along 
the dominant axis (since an ISP located at the opposite high-polarization end of the data 
distribution would be inconsistent with the observed data of the first epoch). The data of 
the second epoch limits the ISP Stokes parameters to Qisp < 0.3% and Uisp < 0.1%. 
The observation at the third epoch shows the SN in the nebular phase, when the density of 
scattering particles is extremely low. At this epoch it is expected that the density is insuf- 
ficient to produce a net intrinsic polarization, rather that the observed polarization at this 
epoch is just due to the combined ISP of the foreground and host galaxy ISP. The average 
values of Q and U data would therefore be equivalent to Qisp and Ujsp- In order to use 
this data to determine an ISP, the data were rebinned to 50A to improve the S/N. The 
polarization signature at the red end of the data suffers from low S/N, due to the response 
function of FORSl at those wavelengths. This leads to a polarization signature with larger 
uncertainties. This section of the spectrum, with A > 8OOOA, was removed from the calcu- 
lation. The rebinning of the spectrum, and removal of high- uncertainty wavelength range, 
significantly reduced the scatter in the data and the centroid of the data on the Q — U plane 
was measured, using separate weighted averages for both Q and U, to be Q = 0.12 ± 0.07 
and U = —0.12 ±0.08. The quoted uncertainties on these values are the standard deviations, 
which provide a better estimate of the scatter than the standard error on the mean. These 
values are consistent with the limits placed on them using the observations from the first 
and second epochs, and are the same as the measured Stokes parameters in the regions of 
the spectrum of the t hird e poch with the highest S/N (see §3.2.11) . 

The ISerkowski et al.l (1l975l ) form of the ISP has only a small degree of A-dependence, across 
the wavelength range being studied here. The small value of the determined Stokes pa- 
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rameters of the ISP are sufficiently low that any effects due to the A-dependence of the 
ISP would be much less than the quoted uncertainties. We adopt, therefore, values of the 
ISP Stokes parameters of Qisp = 0.12 ± 0.07 and Ujsp = —0.12 ± 0.08, corresponding to 
Pisp = 0.17 ± 0.08 and 9jsp = 157?5 (^ yhich is in the direct ion tangential to the host spiral 
arm in NGC 7424, as shown on Fig. [U IScarrott et al.lll987l ). These values limit the redden- 
ing, assuming a Galactic polarization law, to E{B — V) > 0.01 as a lower limit, requiring a 
reddening component in addition to the Galactic foreground reddening. 
Correction of the ISP requires the vectorial subtraction of Q and U, effectively an offset of 
the zero-point of the intrinsic Q and U parameters from the observed Stokes parameters. The 
principal ramification of this property is that the corrected total polarization po ^ Pobs—PisPi 
and polarized features can become depolarized and vice versa. The polarization spectra, af- 
ter correction for the ISP, are also shown for the three observational epochs as Figs. [31 IH 
and El 



Howell et al.l (120011 ) suggest that the region between 4800-5600A should be depolarizing. 



since in that wavelength range the line blanketing opacity of Fe II lines dominates over 
electron scattering, leading to net depolarization. We note that in the case of SN 2001ig, 
however, the observed polarization over much of that wavelength range is similar in magni- 
tude to that observed for Ha for which the presence of net intrinsic polarization has already 
been discussed. The narrower wavelength region of 4600-4800A shows lower levels of polar- 
ization, at the ffist epoch, than 4800-5600A. Only in the narrower region 4600-4800A is there 
almost complete depolarization evident, with the observed polarization consistent with the 
ISP alone (see Fig. [6]). The Fe II lines are further discussed in §3.2.61 



3.24. Ha, St II and He I 6678 A 

At the ffist epoch, the Ha polarization signature seems counter-intuitive, with the max- 
imum polarization of this feature observed at 6350A, near the emission peak in the flux 
spectrum, and minimum polarization reached just blueward of the absorption minimum. 
Low polarization is generally seen for the entire absorption feature (see Fig. [3|), with the 
Stokes parameters lying on or close to the dominant axis with no obvious loop structure 
(Figs. [9] and [TTal) . At day 31, the blue edge of the Ha/Si II absorption proffie (~ 6IOOA) is 
coincident with a peak in polarization of 0.7% (Figs. HI and [TU]) . Significant polarization is 
observed at both of the flux minima of Ha and Si II, but it is smaller in magnitude than the 
polarization at the blue edge. Redward of the flux absorption minima is a partly depolarized 
region at 6300A associated with the Ha emission. The polarization has a local minimum 
of 0.2% at the location of the He I 6678A emission component. Within this depolarized 



- 12 - 



wavelength range there is, however, a peak in polarization at ~ 6400A, comparable to that 
observed for the blue edge of the Ha/Si II absorption profile. This peak coincides with the 
absorption minimum of the He I 6678A P Cygni profile, the prominent "notch" at 6400A 
in the fiux spectrum (Figs. HI and [TU]) . This polarization signature across the He I 6678A 
absorption minimum matches the polarization peak observed for the absorption minimum of 
He I 5876A, where the peak polarizations are equal and the Stokes Parameters are observed 
to rotate through 10°. 

The behaviour of the blended Ha/He I/Si II feature on the Q-U plane, presented in Fig. 
I12al shows that data at these wavelengths are aligned with the dominant axis, but that the 
H I absorption feature forms a loop deviating from a straight line by ±0.2%, in the direction 
orthogonal to the dominant axis. The rotated Stokes parameters at this epoch (Fig. [TOl) 
actually resolve the separate structures of the blend. A peak in the orthogonal rotated Stokes 
parameter is associated with the absorption minimum of Ha in the fiux spectrum, but there 
is no corresponding feature in P^- The absorption feature associated with He I 6678A giving 
rise to the "notch" in the blended profile in the fiux spectrum, is, conversely, only observed 
in Pd and not in Pq. The loop structure observed on the Q — U plane for this blend is, 
therefore, due to the superposition of two lines arising from different species with different 
geometries within the ejecta. By day 256, Ha nebular emission is not significantly detected 
in the fiux spectrum and there is no significant polarization at the expected wavelength for 
this line. 

3.2.5. I 7774 A 

The O I 7774A feature is observed in the fiux spectra at all epochs, and in the case of 
the last epoch it is only seen in emission. At the first epoch there is no excess polarization, 
over the wavelength range of the O I 7774A line (which is observed at the second epoch) 
as shown in Fig. Illbl At 2002 Jan 3, the standard polarization profile expected for P 
Cygni spectroscopic profiles is observed, with peak polarization (0.7%) associated with the 
absorption component at ~ 7500A and decreasing levels of polarization to a minimum of 
0.3% at the peak of the emission component (Fig. Hj). These maxima and minima are 
similar to the observed polarization levels of Ha and He I 6678A. The presence of a telluric 
absorption band at the wavelengths of this feature complicates both the spectroscopic and 
spectropolarimetric interpretation of this line. There is sufficient resolution, even with the 
data rebinned, to safely identify those bins with high uncertainties on the Stokes parameters 
and exclude them from the analysis. On the Q — U plane, on Fig. I12bl the O I data points 
are tightly clustered about the dominant axis, with a maximum deviation in the orthogonal 
direction of ±0.17%. The O I 7774A is almost completely described by P^, with nearly no 
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signal or variability in Pq over the entire line profile (Fig. [T0|) . At day 256, 2002 Aug 16, 
the line is observed in emission, and there is no significant polarization signature associated 
with it. 



3.2.6. Fe II 



The wavelength range from 4800-5400A is expected to be mostly composed of Fe II 
lines, with their o verlap (due to Dop pler broadening) leading to depolarization across this 
wavelength range (IHowell et al.ll200lh . At day 13, there is a net polarization of 0.2% observed 
across this range. This is likely due to the presence of strong polarization signatures from 
other lines, due to other species, which cannot be completely depolarized by the Fe II lines; 
for example, Fe II 4924A blended with H/? and He 4921A, and 5018A blended with He 5015A 
( iMatheson et al.ll2000al ). The Stokes parameters show no preferred orientation across the Fe 
II lines (Figs. IH] and lllcp and a rise in the dominant rotated Stokes parameter across the 
wavelength range is interpreted as being due to blending with other species. At day 31, a 
net depolarization from 0.5% to 0.1%, compared to other features, is observed across this 
wavelength range. The effects of possible blending are illustrated by the proximity of H/5, 
blueward, and He I 5876A redward. Taking the wavelength range of 4800-5400A as being 
representative of a number of Fe II lines, the data on the Q — U plane and the dominant 
and orthogonal rotated Stokes parameters shows that the associated polarization is well 
described by the dominant axis (see Figs. [10] and I12cl) with minor deviations due to the 
behavio ur specific to individua l lines of either Fe II or blends of Fe II lines with those other 
species ( IMatheson et al.ll2000al ). The cause of the polarization peak at the second epoch at 
4750A is unclear, given the high density of Fe II lines and H/3. The polarization feature 
is significantly broader than any of the individual absorption features, in that wavelength 
range, observed in the flux spectrum. This implies that it is unlikely that a single line is 
producing it. This polarization is only observed along the dominant axis (Fig. [T0]and ri2cp . 
unlike Ha. It is, therefore, unlikely to be due solely to H/3 and since it would be expected to 
have similar polarization properties as Ha. It is more likely to have arisen from Fe II lines 
as it is oriented along the dominant axis, similar to the polarization properties measured for 
the Fe II lines at larger wavelengths. 



3.2.7. Ca II IR triplet 



The polarization properties of the Ca II IR triplet absorption feature suffer from a 
number of effects which increase the levels of uncertainties on the measurements. The region 
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redward of 8000 A shows a significant drop in the response of the FORSl detector and hosts 
a particularly broad telluric absorption band. At the first epoch, there is a polarization 
associated with the absorption feature, that reaches a maximum degree of polarization of 
0.4 ± 0.3% at the absorption minimum (see Fig. Illdl) . At day 31, with increased S/N 
in the red, due to the increased brightness of the SN, the maximum polarization is again 
observed to directly correlate with the absorption minimum and is measured as 0.7 ± 0.2%. 
As for the O I 7774A line, the uncertainties of the measured Stokes parameters are large 
and this is due to the telluric absorption band. At the red extreme, the profile of the Ca 
II IR blended triplet absorption, in the fiux spectrum, appears rising towards the emission 
component. The polarization angle is observed to rotate from the continuum polarization 
angle of ~ 40° through +20° moving redward towards the absorption minimum and back 
through —20° when moving from the absorption minimum toward the unobserved emission 
component. On the Q — U plane, see Fig. I12dl an incomplete loop is observed, as the emission 
component was not observed. The loop appears to be aligned with the direction orthogonal 
to the dominant axis (deviating by a maximum 0.4%), and is approximately parallel to the 
dominant axis determined for the observations of day 13. The rotated Stokes parameters Pd 
and Po show that at day 13 the polarization of the Ca II is parallel to the dominant axis 
(Fig. [9l). By day 31, however, the polarization has substantial components in both P^ and 
Po (Fig. [TOl) . At day 256, there is no significant polarization associated with the Ca II IR 
triplet, with large error bars being symptomatic of the low S/N at these wavelengths. 



4. Discussion 



The observations of SN 200 lig show the three distinct epochs in the spectroscopic evo- 
lution of a Type lib SN, accompanied with three distinct sets of polarization properties. 
The fiux spectra of SN 2001ig at all three epochs are compared on Fig. [2] with spectroscopic 
observations of SN 1993J, obtained from the SUSPECT archiv^. At the earliest epoch, 
SN 2001ig is observed to have a higher ratio of Balmer to helium line strengths than SN 
1993J. Indeed, at this epoch the sp ectrum of SN 2001ig is more similar to that of SN 1987K 
( Filippenkol Il988l : IClocchiattil 12002 ) . with the absorption of He I 6678A at higher velocities 
and, therefore, being observed to be blueward of Ha. The absorption minimum of Ha how- 
ever, is observed to be significantly more blueshifted, by ~ 3500km s~^ , than observed for 
1987K. Similarly, at 31 days the absorption of He I 6678A is at hi gher velocities than was 
observed for SN 1993J at a similar epoch (12000 vs. 6940km s~^ ; ISwartz et al.l 119931 ). In 
the nebular phase, at day 256, the spectrum of SN 2001ig is very similar to SN 1993J at 



■"^http: / /bruford. nhn.ou.edu/~suspect/ 
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a similar age showing [Ca II] 7291,7324A, He I 5876A and O I 6300/6364A. The ratios of 

ines in SN 2001ig are approximately similar to those obser ved for SN 



the strengths of these 
1987K and SN 1993J ( 



Fi 



ipDenkolll988l : iMatheson et al.ll2000al : iRvder et al.ll2006f ). although 



Filippenko fc Chornockl fl2002h note that the Mg I] 447lA observed for SN 2001ig is partic- 



ularly strong (see Fig. [2]). 

At 13 days, when hydrogen features dominate the spectrum, the general spectrum-wide po- 
lariza tion is low ~ 0.2%. This is consistent w ith the low polarizations observed for Type IIP 
SNe (IWang et al.l l2002l : [Leonard et al.l l200ll ) which have extended optically-thick hydrogen 
envelopes. The presence of a non-zero polarization implies that there is an asymmetry of 
the outer hydrogen layers at the early times * which for an edge-on oblate spheroidal config- 
uration only requires an axis ratio of ~ 0.9 (jHoflichlll99ll ). 

It is only at the second epoch when the expected depolarized Ha emission feature is first 
observed. This is due to the decrease in the velocities of the ejecta material compared to 
the first epoch, which changes the way in which the absorption and emission features of Ha 
and He I 6678A are blended. The increase in polarization of SN 2001ig between the first 
and second epochs is also similar to other Type II SNe, when the hydrogen layers c ease to 
be op tically thic k and the interior as ymmetric He core is revealed (e.g. SN 1987A, iJeffery 
199ll : SN 2004di. [Leonard et al.ll2006[ ). 

The evolution in the polarization signature of SN 2001ig at the first and second epochs is 
matched by the spectroscopic evolution of the SN, with the increase in the strength of helium 
lines by the second epoch. The data on the Q — U plane, at this epoch, are aligned with 
an obvious dominant axis which suggests that the He core layer asymmetries still re tain a 
strong axial symmetry although clearly not a spherical symmetry (jWang et al.ll2003bl ). The 
dominant axis determined for the data at the second epoch is, however, almost orthogonal 
to the dominant axis of the data at the first epoch. This suggests that while asymmetries 
were present at day 13, they were not coupled to the larger asymmetries present in the core 
that were revealed by day 31 in either the degree of the asymmetry or the orientation of 
the asymmetry on the sky. This provides a case, therefore, that the hydrogen envelope was 
sufficiently extended to not be di rectly coupled to the asymmetries of the He core layers 
due to the explosion mechanism. iHofiichI (119951 ) suggested that the small degree of early 
polarization observed for SN 1993J was due to tidal distortion of the progenitor by a binary 
companion, if the SN arose in such a system. Alternatively, or in addition, another source 
of polarization at the earlier epoch might be the aspherical distribution of sources of exci- 
tation, from newly synthesized material, causing the shape of the photosphere to deviate 
from a spherical symmetry, while the density distrib ution of the hydrogen envelope remains 
physically spherically symmetric ( Hoflich et al.ll200ll ). The asymmetry measured at the first 
epoch is, perhaps, an upper limit on the physical asymmetry of the hydrogen envelope. 
SN 2001ig and SN 1993J exhibit the same evolution in polarization properties between days 
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13 and 31, with an increase in the degree of the polarization accompanied by the rotation 
of the polarization axis. At the first epoch, the l^-band polarization of SN 200 lig (from 
synthetic broad-band polarimetry) is 0.3% with 6 = 185° (ISP subtracted), whereas SN 
1993J was polar ized at a level of ~ 0.6% at approximately the same epoch (1993 Apr 7; 
Tran et al.lll997l ). At 31 days, the broad-band V pola rization for SN 2 001ig is 0.6%, whereas 



the polarization of SN 1993 J was 0.8% (1993 Apr 26; iTran et al.lll997l ). and the polarization 
angle had rotated through ~ 40°. In absorption components, however, SN 2001ig is less 
polarized than SN 1993J by <^ 1%. The meandering p ath of the V-ban d polarization of SN 
1993 J, as plotted on the Q — U plane (see Fig. 10 of iTran et al.l 119971 ) . is also seen for SN 
2001ig; although the number of observational epochs of SN 2001ig is much lower. 
In these respects SN 2001ig is similar to 1993J, but there are important exceptions. At 
day 31 the polarizations of the He I 5876A and 6678A absorption minima of SN 2001ig are 
comparable to that of the Ha absorption feature, whereas there is no observed polarization 
signature attributable to He 6678A in SN 1993J. This is due to the He I 6678A line for SN 
200 lig being both stronger and having a higher velocity than was observed for the same 
line in SN 1993J. The higher velocity for the absorption of He I 6678A means that it is not 
completely coincident with the depolarizing Ha emission flux, such that the polarization 
signature of the ab sorption is still observed. 

Wang et al.l (l200ll ) suggested that the classification of a CCSN as being of Type lib might 
depend on the orientation of the SN as seen by the observer, given the spectropolarimetric 
similarity between SN 1993J and SN 1996cb. The observations presented here have shown 
that SN 2001ig arises from the same family of SNe as SN 1993J and SN 1996cb and, by 
extension due to spectroscopi c similarity, to SN 1987K, but that Type lib SNe are not as 
homogeneous as suggested by lWang et al.l (120011 ). In the context of Type I CCSNe, we note 
that a whol e range of polarizations have been measured for Type Ib/c SN e, from ~ 1% for 
SN 2002ap jMazzah et al.l boosi : IWang et al.l l2003bl : lleonard et aP l2002h to ^ 4% for the 
Type Ic SN 1997X (although this is a broad band polarization measurement and it is un- 
clear if this value is representative of the continuum polarization; IWang et al.ll200ll ). On the 
other hand Type IIP and lib SNe exhibit a more limited range in polarization than the more 
hydrogen deficient SNe. This suggests that the retained hydrogen envelope of the progenitor 
is important to dictating the polarization evolution of Type lib SNe, and that Type lib 
SNe which have undergone the transition to the hydrogen poor phase are not necessarily the 
same as "pure" Type lb SNe. 

Particular note is made of the loop structure observed, on the Q — U pl ane, for certain spec- 
tral lines. Loops were observed in spectropolarimetry of SN 1987A (ICropper et al.l Il988l ) 
and similarly have been presented for the Type la SN 2001el (IKasen et al.l 120031 ). The pro- 
duction of a loop requires an increase (or a decrease) in polarization across a spectral line, 
accompanied by a rotation in polarization angle. While IWang et al.l (l2003bf ) discussed the 
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"dominant" and "orthogonal" loop is a mixture of both. The rotated Stokes pa- 

rameters of the Ca II IR triplet cannot be decomposed at day 31 to two distinct orthogonal 
components (Fig. (TU]), leading to loop structure on the Q — U plane (Fig. I12al) . The Ha/He I 
blend is not, however, a true loop since it can be decomposed into two orthogonal components 



depen ding on which element dominates absorption at a particular wavelength. iKasen et al. 



(120031 ) discussed how loops may be formed by breaking a global axial symmetry by having 
varying degrees of the departure of the ejecta from an axial symmetry as a function of depth 
or with a clumpy SN ejecta. The loops may hint at the presence of plurne-like structure in 
SN 2001ig e j ecta, which is observed in SN remnants (IFesen fc Gundersonlll996l ). 



Wang et al.l (120061 ) applied straight line fits to the Stokes parameters on the Q — U plane to 



judge the relative distributions of particular elements within the ejecta of the Type la SN 
2004dt. If spectral lines produce loops rather than a single dominant axis (i.e. a straight 
line) on the Q — U plane, then straight-line fits of a dominan t axis to a lo op would be in- 
appropriate; this may then explain the high values of xl that IWang et al.l found for fitting 
dominant axes to a number of spectral features of SN 2004dt. 

The locations of loops, hence the range of polarization angles, on the Q — U plane can be 
used in a similar manner to fits of the dominant axes to particular spectral features to look 
for the relative distributions of different elements within the SN ejecta. For day 31, species 
such as He and O are observed to have similar orientations on the Q — U plane and all lie 
along the dominant axis (see Figs. [TO] and [T2cl) . suggesting a similar distribution for these el- 
ements in the ejecta. The polarization of Ha, however, is predominantly orthogonal to these 
species and, hence, has a significantly different distribution within the ejecta. The Stokes 
parameters across the Ca II absorption feature, however, seem to have a similar orientation 
to the dominant axis measured for day 13. 

The nature of the proge nitor system, single or binary, that gave rise to SN 2001ig has been 
subject to much debate. iRyder et al.l (120041 ) observed perio dic fluctuatio ns in the radio light 
curve of SN 2001ig, due to fluctuations in the CSM density. iRyder et al.l concluded that this 
was due to the presence of a binary companion, disturbing the stellar wind from the pro- 
genitor as the immediate CSM was formed, with the period i city i n the density fluctuations 
corresponding to the orbit of the companion. iKotak fc VinkI (120061 ) suggested, however, that 
episodic mass loss from a Luminous Blue Variable phas e of the progenito r could also produce 
similar density fluctuations in the CSM. Additionally ISoderberg et al.l (120061 ) showed that 
SN 2003bg had a similar radio light curve to SN 2001ig; the likelihood of two systems having 
the same orientation and binary progenitor system parameters, to produce such similar light 
curves, was considered to be too low to make a binary progenitor a plausible scenario based 
on radio measurements alone. In the case o f the Type lib SN 1 993 J, the progenitor was 
identified in fortuitous pre-explosion images (lAldering et al.lll994) and the companion was 



identified in very late time images and spectroscopy (iMaund et al. 



20041 ). demonstrating the 
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produ ction of that particular Type lib SN from a binary progenitor system. iRyder et al. 



( I2OO6I ) have claimed to observe a possible stellar source, 3 years post-explosion, at the site of 
the fading remnant of SN 200 lig. It is, therefore, particularly interesting to dis cuss the ob 



servations presented here in the context of the nature of the progenitor system. IWang et al. 



( I2OOII ) concluded that Type lib SNe, having so many similar spectropolarimetric properties, 
must arise with similar orientations to the observer. The shared characteristics of SN 200 lig 
with Type lib SNe 1993J and 1996cb, while not a complete match, might suggest a common 
type of progenitor with similar orientation; thus favouring a binary progenitor. The differ- 
ences in the polarization properties observed for SN 2001ig days 13 and 31 suggests that the 
hydrogen and core layers did not possess the same asymmetry, both in degree and orienta- 
tion. The hydrogen envelope was, therefore, sufficiently large (in radius) t o not be subject 



to the same asymmetries as the core; this is in contrast to the picture of ISoderberg et al. 
( l2006l ). who suggest a thin hy drogen layer on a W olf-Rayet progenitor. 
In both the je t-torus model (IHoflich et al.l I2OOII ) and the binary progenitor system model 
( iHoflichl I1995I ). it is clear how different asymmetries of the envelope and the core- layers 
(defined by the physically orthogonal equatorial plane and rotational axis of the progeni- 
tor) would lead to different polarization angles for the dominant axis between the first two 
epochs. In the jet-torus model, the production of Ni in a jet would cause asymmetric exci- 
tation of the hydrogen envelope aligned with the rotational axis, while the core-layers would 
form a torus in the equatorial plane. In the binary progenitor model, the asymmetries of 
the hydrogen envelope would be caused by tidal interaction with the companion, causing 
deformation in the direction of the equatorial/orbital plane, and a jet, that did not breach 
the outer core-layers, would cause elongation of the core-layers along the rotational axis. 
Neither of these two simple models can explain the physical rotation of the asymmetries, 
between the fi rst and second epoc hs, of only ~ 40° when both favor a rotation of 90°. 
Alternatively, iMaund et al.l (120071 ) present a "tilted jet" model. In this model the jet axis 
is not aligned with the rotational axis of the progenitor, such that the outer layers and the 
inner core layers do not have to have either the same or orthogonal axes of symmetry. The 
outer layers, in this model, retain the shape of the progenitor, but the core-layers are elon- 
gated by the mechanical force of the jet and the redistribution of Ni along the jet axis. In 
the tilted-jet model, the angle between the axis of symmetry of the core and the outer layers 
is not fixed. We note that the pola rization axis of SN 1987A is tilted by ~ 15° with respect 
to that of the circumstellar rings (jWang et al.l |2002| ) and similar deviations might charac- 
terize Cas A ( IWheeler et al.l 120071 ). The Standing Accretion Shock Instability in 3-D may 
provide a mechanism through which the axis at which the jet is launche d is no longer aligned 
with the rotational axis of the progenitor, with lBlondin fc Mezzacappal (120071 ) modelling the 
rotational axis of the proto-neutron star, in 3-D core-collapse simulations, becoming mis- 
aligned by angles such as 10°, 15° and 45° from the rotational axis of the Fe core. Other 
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non-axisymmetric instabilities may also play a role (jWheeler fc Akiyamal 120071 ). 
It is clear, however, that along with appropriate models spectropolarimetry data can provide 
important insight into the nature of the explosion mechanism and the progenitors of this 
type of CCSNe. 



5. Conclusions 



The Type lib SN 2001ig has been observed to have a strong polarization signature, 
consistent with other Type lib SNe. The polarization properties are generally in excess of 
those observed for Type IIP SNe, but at the lower limit of polarizations observed for the 
hydrogen deficient Type Ibc SNe. The ISP is constrained by analyzing the polarization 
behaviour of SN 2001ig at three separate epochs, but particularly using observations at the 
last epoch, in the nebular phase, when SN 200 lig is believed to be intrinsically unpolarized. 
SN 2001ig has a low continuum polarization of 0.3% at 13 days, consistent with an almost 
spherical photosphere (with deviation of ^ 10%) and the spectrum is dominated by hydrogen 
features. At 31 days the decrease in the strength of the hydrogen lines and the increase in 
the strength of the helium lines, along with a sharp increase in continuum polarization to 
~ 1%, demonstrate that at this epoch the highly asymmetric He core is being revealed. 
The absorption components of He I 5876A and 6678A are observed at 31 days to be highly 
polarized (to ~ 0.8%). Loops, on the Q — U plane, are observed for Ha/He I 6678A, O I 
7774A and the Ca II IR triplet at 31 days. The presence of these loop features indicates a 
deviation of the SN ejecta from a global axial symmetry. There is significant rotation of the 
observed Stokes Q and U parameters, between the first and second epochs, through an angle 
of ~ 40°; indicating that the hydrogen envelope was sufficiently extended to not be coupled 
to the asymmetries of the He core-layers, associated with the explosion mechanism. The 
different geometries of the outer envelope layers and the inner core region may be evidence 
of the "tilted jet" model. SN 2001ig is intrinsically unpolarized at 256 days. The different 
polarizations at the three epochs show the same physical process observed in Type IIP SNe, 
of an almost spherically symmetric hydrogen layer shielding a highly asymmetric He core. 
As time since explosion increases the hydrogen layer becomes optically thin revealing the He 
core and the polarization is seen to increase. The expansion velocities of the hydrogen and 
helium in SN 2001ig are significantly higher than observed for 93J-like Type lib SNe, but are 
similar to those observed for SN 1987K, and the polarization properties of SN 2001ig show 
distinct differences from SN 1993J and SN 1996cb, which are mostly due different degrees 
of line blending. The differences between SN 2001ig and SNe 1993J and 1996cb suggest 
that Type lib SNe are not as homogeneous as perhaps previously considered, but may be 
sufficiently homogeneous to share similar geometries and, potentially, progenitor systems. 
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Table 1: Journal of Spectropolarimetric Observations of SN 2001ig 



Star 


Date 


Exposure 


Wavelength 


Air Mass 


Filter 


Gomment 






(s) 


(A) 








SN 2001ig 


2001 Dec 16.11 


4 X 1000 


4450-8635 


1.86 


GG435 


SN 


HD 49798 


2001 Dec 16.12 


10 


4450-8635 


1.31 


GG435 


Flux STD 


LTT 9491 


2002 Jan 03.02 


200 


4450-8635 


1.47 


GG435 


Flux STD 


SN 2001ig 


2002 Jan 03.07 


4 X 500 


4450-8635 


2.19 


GG435 


SN 


EG 274 


2002 Aug 15.98 


60 


3700-8635 


1.03 


None 


Flux STD 


SN 2001ig 


2002 Aug 16.26 


4 X 1200 


3700-8635 


1.04 


None 


SN 


GD 50 


2002 Aug 16.44 


60 


3700-8635 


1.09 


None 


Flux STD 



All observations conducted with the GRIS-300V grism. 



Table 2: Prominent P Gygni Profiles and Velocities at Absorption Minimum at -|-13 and -|-31 

days 



Epoch 
-1-13 days -1-31 days 

Line 





V (km s ) 


V (km s ^ ) 


Fe II 5169A 


blended 


-10 300 


He I 5876A 


-14 000 


-12 000 


Si II 6355A 




-8 000 


Ha 


-19 000 


-13 600 


He I 6678A 


-19 000 


-12 000 


I 7774A 


-12 700 


-10 600 


Ga II IR triplet 


-12 700 


-12 400 
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Fig. 1. — Digital Sky Survey image of the location of SN 2001ig, shown by the filled black 
circle and cross hairs, relative to its host galaxy NGC 7424. The arrow shows the orientation 
of the Interstellar Polarization component determined in Section I3.2.3[ 
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Fig. 2.— Flux spectra of SN 2001ig at 2001 Dec 16, 2002 Jan 3 and 2002 Aug 16 (or 13 
31 and 256 days post-explosion respectively). Line identifications follow iMatheson et al. 
(j2000al ). Telluric lines are indicated by the ® symbol. Also shown are spectra of SN 1993 J, 
from the SUSPECT archive, similarly scaled for sim ilar post-explosion epochs, as given on 
the figure, using the explosion date of SN 1993 J of iTran et al.l (119971 ). The vertical lines 
show the rest wavelengths and, where appropriate, the associated absorption profile. The 
wavelengths of these spectra have been corrected for the recessional velocities of the host 
galaxies of these two SNe. 
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Fig. 3. — Spectropolarimetry of SN 2001ig, acquired on 2001 Dec 16 (13 days post-explosion), 
showing the total observed polarization p, the total polarization after correction for the ISP 
Pq, the Stokes Q and U parameters, the polarization angle 9 and the total flux spectrum 
(ergss~^ cm~^ A~^). The Stokes parameters have been re-binned to IsA. The thick grey 
lines indicate the position of the ISP Stokes Q and U parameters, effectively the zero-point 
for the Stokes parameters intrinsic the SN. The zeropoint for observed Q and U is indicated 
by the horizontal dotted line. All data is presented with the wavelength scale corrected for 
the heliocentric recessional velocity of the host galaxy. 
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'2002-01-03 




5000 6000 7000 8000 

Rest Wavelength / A 

Fig. 4. — Same as Fig. [3l but for spectropolarimetry of SN 2001ig acquired on 2002 Jan 3 
(31 days post-explosion). The data is uncorrected for the ISP. Note the small dispersion of 
the polarization angle 9, compared to the earlier epoch of Fig. [3l 
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Fig. 5. — Same as Fig. [3l but for spectropolarimetry of SN 2001ig acquired on 2002 Aug 16 
(256 days post-explosion). 
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Q(%) 

Fig. 6. — Stokes Q and U parameters, as a function of wavelength from observations acquired 
on 2001 Dec 16 (13 days post-explosion). The Stokes parameters have been re-binned to 
50A. The data have been corrected for the ISP (shown as the solid black circle; see text). 
The dominant axis (see text) is indicated by the straight line. The wavelength of each point 
is indicated by the colour, following the scheme of the colour bar on the right hand side. 
The location of the origin of the Stokes plane {Q = and U = 0) is indicated by cross. 
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Fig. 7. — Same as Fig. [6], but for spectropolarimetry of SN 2001ig acquired on 2002 Jan 3 (31 
days post-explosion). Note that the dominant axis has rotated physically ~ 40° compared 
to that of the earher epoch of Fig. El 
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Fig. 8. — Same as Fig. El but for spectropolarimetry of SN 2001ig acquired on 2002 Aug 16 
(256 days post-explosion). 
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Fig. 9. — The Stokes parameters, corrected for the ISP, rotated to the dominant (P^) and 
orthogonal axes (Pq) for the data of 2001 Dec 16, compared with the flux spectrum at this 
epoch. There is little difference between these rotated Stokes parameters and the Q and U 
Stokes parameters (Fig. [3]), since the dominant axis is nearly aligned [6^ = 4.5°) with the Q 
axis. Zero polarization is indicated by the heavy grey line. 
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Fig. 10. — The Stokes parameters, corrected for the ISP, rotated to the dominant {Pd) and 
orthogonal axes (Po) for the data of 2002 Jan 03, compared with the flux spectrum at this 
epoch. Zero polarization is indicated by the heavy grey line. 
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Fig. 11a. — Loops on the Q-U plane for Ha/He I 6678A blend, the wavelength range covering 
O I 7774A, Fe II lines in the range 4800-5600A, and Ca II IR absorption from the observations 
of 2001 Dec 16. The data have been re-binned to 15A. The heavy dashed line indicates the 
dominant axis, calculated for the entire data set for this epoch shown as shown on Fig. El 
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Fig. lid.— 
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Fig. 12a. — Loops on the Q-U plane for a) Ha/He I 6678A blend (Points of interest are 
indicated by the numbers, following the wavelength progression of the data on the Stokes 
plane: 1) Blue edge of the Ha absorption profile in the flux spectrum; 2) Absorption minimum 
of Ha, 3) Absorption minimum of He I 6678A; and 4) The Ha and He I 6678A depolarizing 
emission feature), b) O I 7774A, c) Fe H hues in the range 4800-5600A, and d) Ca II IR 
absorption from the observations of 2002 Jan 3. The data have been re-binned to 15 A. The 
heavy dashed line indicates the dominant axis, calculated for the entire data set for this 
epoch shown as shown on Fig. [71 For Ca II the dominant axis determined at the first epoch 
is shown as the dot-dashed line. 
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Fig. 12c.— 
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